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Abstract 

The influence of three well-known disaccharides, namely trehalose, maltose and 
sucrose, on some structural and dynamical properties of lysozyme has been inves- 
tigated by means of molecular dynamics computer simulations in the 37-60 wt % 
concentration range. The effects of sugars on the protein conformation are found 
relatively weak, in agreement with the preferential hydration of lysozyme. Con- 
versely, sugars seem to increase significantly the relaxation times of the protein. 
These effects are shown to be correlated to the fractional solvent accessibilities of 
lysozyme residues and further support the slaving of protein dynamics. Moreover, a 
significant increase in the relaxation times of lysozyme, sugars and water molecules 
is observed within the studied concentration range and may result from the perco- 
lation of the hydrogen-bond network of sugar molecules. This percolation appears 
to be of primary importance to explain the influence of sugars on the dynamical 
properties of lysozyme and water. 

Keywords : Biopreservation, disaccharides, preferential hydration, protein stability, hy- 
drogen bonds 



2 



1 Introduction 



The preservation of biological molecules like proteins is a fundamental goal of pharmaceu- 
tical, food, and cosmetic industries [1-4]. Indeed, they generally remain stable only under 
very stringent conditions of temperature, pH, hydration level or ionic strength [3,5]. The 
addition of stabilizing agents is a common strategy to increase the stability of proteins 
both in solution or in the dry state [3,5]. Among these compatible solutes, trehalose (see 
figure la), a disaccharide (C^-f^On), has been found particularly efficient [2,6]. This 
sugar is naturally synthesized in large amounts by plants, insects and microorganisms 
such as yeasts and nematodes which are able to withstand for an extended period severe 
environmental conditions of low/high temperatures and/or drought experienced in desert 
or polar regions [7,8]. These organisms enter into a biological state where their activity 
is almost completely suspended as long as environmental conditions remain deleterious, 
and then resume their normal activity. This phenomenon is called anhydrobiosis. Despite 
the many experimental [6,9-11] and numerical [12-14] studies which have been devoted 
to the bioprotection phenomenon, the molecular mechanisms at its origin are still not 
clearly understood. Several hypotheses, mostly based on the properties of trehalose - 
high solubility, low reactivity, good glass-forming and antioxydant properties, etc. - have 
been proposed to explain why it is the most effective bioprotectant among sugars and 
polyols [8]. Nonetheless, none of them is fully satisfactory, since it generally covers only 
a limited range of temperatures and hydration levels. 

In aqueous solutions proteins may be thermodynamically stabilized by the preferential 
exclusion of cosolvent molecules from the protein/solvent interface, which make thermo- 
dynamically unfavorable the increase of the solvent accessible surface area (SASA) of 
proteins upon denaturation [19]. This exclusion could arise from excluded volume effects 
and an increase in the surface tension of water induced by solutes [20]. This may explain 
why among many solutes trehalose was found to be the most excluded [21,22], given its 
larger hydration number in comparison to other disaccharides such as maltose and su- 
crose [23,24] (the structures of these sugars are shown in figures lb and c, respectively). 
Moreover, solutes may hinder the formation of ice, which is generally lethal to organisms. 
Magazu et al. have suggested that trehalose would be more effective in inhibiting ice for- 
mation than other disaccharides, because it binds to a larger number of water molecules 
and would therefore have a greater destructuring effect [25]. 

At high temperatures or high osmotic pressures, dehydration is observed. Biolog- 
ical molecules may experience important stresses when hydration water molecules are 
removed [26]. The effect of the addition of stabilizing solutes may be two-fold. First, 
the solution is likely to vitrify before the complete removal of hydration water molecules, 
since the solution viscosity raises exponentially when decreasing the water content. The 
glass formed may kinetically maintain the conformation of proteins and prevent the fusion 
of membranes. Trehalose would then be a better bioprotectant than maltose or sucrose 
given its higher glass transition temperature T g , as suggested by Green and Angell [27]. 
This hypothesis is generally well accepted, even though counter-examples have been re- 
ported [11,28,29]. The second effect of an increase of the concentration of stabilizing 
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Figure 1: Schematic representation of trehalose (a), /3-maltose (b) and sucrose (c) in 
their most stable cristalline forms, which are dihydrate [15], monohydrate [16] and an- 
hydrous [17] crystals, respectively. Oxygen, carbon and hydrogen atoms are displayed in 
black, grey and white, respectively. Intramolecular hydrogen bonds (HBs) are represented 
by dashed cylinders. In these conformations, trehalose, maltose and sucrose form zero, 
one and two intramolecular HBs, respectively. This figure was generated using VMD [18] 
(http:/ /www. ks.uiuc.edu/Research/vmd/). 

solutes upon dehydration is the formation of hydrogen bonds (HBs) between solutes and 
biological molecules. Crowe et al. have proposed that solutes were able to substitute to 
the hydration water of biological molecules, therefore stabilizing their solution structure 
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and functionality [6]. This has been supported by the results from many studies [9-11]. 
There is a lot of debate on the relative importance of the solvent vitrification and the 
formation of solute-biomolecule HBs. Nevertheless, both appear necessary to achieve 
long-term preservation [5,6]. 

Alternatively, Belton and Gil have suggested that the direct interaction of the solute 
glassy matrix with lysozyme is relatively limited, because the hydration water is trapped 
in an amorphous state (non-crystalline ice) by the solute glassy matrix [30]. The role 
of trehalose would be to concentrate residual water molecules close to the protein. This 
would be consistent with the hypothesis of the preferential hydration of proteins in dilute 
to semi-dilute solutions [22]. 

Finally, Cesaro et al. [31] have pointed out the possible role of the interconversion 
between the trehalose dihydrate crystal T 2 h 2 o [15] and the metastable a phase T a [32] 
upon dehydration. This interconversion may make easier the rehydration of biomolecules 
and occurs on time scales compatible with the anhydrobiotic protection, and thus may not 
induce fast changes in the volume or internal pressure of cells, in contrast to ice formation. 

Many simulation studies of ternary systems have been published in the last few 
years [12-14,33-36]. Studies of membrane/sugar /water systems [12,33-36] indicate that 
sugars, and trehalose in particular, interact directly with membranes, and may stabi- 
lize them by forming multiple HBs between several adjacent phospholipids. Furthermore, 
Lins et al. [12] have studied trehalose-lysozyme interactions at a concentration of about 18 
wt % and have shown that trehalose molecules do not expel hydration water. They have 
suggested that the presence of lysozyme induces an aggregation of trehalose molecules, 
which trap a thin layer at the protein surface, in agreement with the suggestion of Belton 
and Gil [30] for solid lysozyme/trehalose/water samples. Similarly, Cottone et al. [37] 
have shown that trehalose is preferentially excluded from the protein surface at 50 and 
89 wt %. They also observed an analogous behavior for sucrose at 89 wt %, which was 
found more excluded than trehalose [13]. Furthermore, Dirama et al. [14] have studied 
lysozyme in a trehalose glass at different temperatures and shown a strong coupling be- 
tween lysozyme and trehalose molecules correlated with the dynamics of protein-sugar 
intermolecular HBs. 

This paper aims at understanding how globular proteins in their flexible native state [38] 
are influenced at room temperature (300 K) by stabilizing solutes at intermediate con- 
centrations, i.e. high enough so that the solvent structure and dynamics are quite 
different from those of pure water, but low enough to allow relatively large-scale mo- 
tions of protein surface atoms. This concentration range is particularly important since 
it may link the protein properties in dilute solutions to those in solid matrices. For 
this purpose, we have probed the influence on some structural and dynamical proper- 
ties of hen egg-white lysozyme, a model globular protein, of three well-known homol- 
ogous disaccharides, namely trehalose [a-D-glucopyranosyl-a-D-glucopyranoside], mal- 
tose [4-0-(o;-D-glucopyranosyl)- / 9-D-glucopyranoside] and sucrose [/5-D-fructofuranosyl- 
a-D-glucopyranoside] (see figure 1), in the 37-60 wt % concentration range by means of 
molecular dynamics computer simulations. This specific concentration range was chosen 
because it corresponds to the range of concentrations where the HB network of the present 
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sugars was shown to percolate in our previous numerical study of disaccharide/ water so- 
lutions [39]. 

2 Simulation details 

Molecular Dynamics simulations of hen egg-white lysozyme (hereafter referred to as 
lysozyme) in sugar/water solutions have been performed using the CHARMM program [40], 
version 29bl. The all-atom CHARMM22 force field [41] has been used to model the pro- 
tein. The CSFF carbohydrate force field [42] has been considered for disaccharides and 
water molecules were represented by the SPC/E model [43]. The production simulations 
were performed in the isochoric-isothermal (N,V,T) ensemble. The length of all covalent 
bonds involving an hydrogen atom as well as the geometry of water molecules were kept 
fixed using the SHAKE algorithm [44], with a relative tolerance of 1CT 5 . A 2-fs timestep 
has been used to integrate the equations of motion with the verlet leapfrog algorithm [45]. 

During the different stages of the simulations, the temperatures have been maintained 
constant with weak coupling to a heat bath (Berendsen thermostat [46]) with a relaxation 
time of 0.2 ps. A cutoff radius of 10 A has been used to account for van der Waals 
interactions, which were switched to zero between 8 and 10 A . A Lennard- Jones potential 
has been employed to represent van der Waals interactions and Lorentz-Berthelot mixing- 
rules have been used for cross-interaction terms. Electrostatic interactions have been 
handled by the particle mesh Ewald (PME) [47] method with k = 0.32 A _1 and the fast- 
Fourier grid densities set to ~ 1/A (48 and 64 grid points in the X/Y and Z directions, 
respectively) . 

The starting structure of lysozyme was obtained from the X-ray crystal structure 
solved at 1.33 A (193L entry of the Brookhaven Protein Data Bank) [48]. Most probable 
charge states at pH 7 were chosen for ionizable residues. The total charge of lysozyme 
(+8 e) was then neutralized by uniformly rescaling the charge of protein atoms, similarly 
to ref . [49] . The disaccharide initial conformations have been deduced from neutron and 
X-ray studies of trehalose [15], maltose [16] and sucrose [17]. The sugar concentrations 
on a protein-free basis are 37, 50 and 60 wt %. These concentrations have been purpose- 
fully chosen based on our previous study of sugar/water solutions [50,51]. Indeed, we 
showed that the relative effect of sugars on water may be distinguished above a thresh- 
old concentration of about 40 wt%. Therefore, their relative influence on lysozyme at 
ambient temperature may be characterized above this concentration. Lysozyme and its 
142 crystallographic hydration water molecules were first placed in an orthorhombic box 
with cell parameters a = b = 46.7 A and c = 62.2 A. Then, disaccharide molecules were 
located and oriented randomly around lysozyme, with minimum sugar-protein and sugar- 
sugar distance criteria, which ensure an isotropic distribution of sugars around lysozyme. 
Finally, water molecules non-overlapping with either lysozyme or sugars were randomly 
added in the simulation box. Initial configurations were minimized in three steps, keeping 
first lysozyme and sugars fixed, then keeping only lysozyme fixed and finally keeping free 
all molecules. This minimized configuration was heated to 473 K in the canonic ensemble 
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during 1 to 3 ns, while maintaining fixed the conformation of lysozyme to prevent con- 
formational changes. This aimed at equilibrating solvent configurations, particularly the 
position and orientation of sugars. Then, the resulting configurations were thermalized at 
300K and simulated in the isobaric-isothermal (N,P,T) ensemble. The stabilized volume of 
the simulation box during this simulation was considered to compute the averaged density 
of the system and used to perform the subsequent simulations in the (N,V,T) ensemble. 
A steepest-descent minimization procedure of 1000 iterations was then performed, whilst 
applying a decreasing harmonic potential on atoms of lysozyme. After the minimization 
procedure, the temperature was raised from to 300 K, with a 5-K temperature increase 
every 250 steps. Then, an equilibration at 300 K was performed during about 80 ps. 
Finally, simulations of 10, 12 and 17 ns were performed for the sytems at concentrations 
of 37, 50 and 60 wt %, respectively, and configurations were saved every 0.25 ps. A 
control simulation of lysozyme in pure water was done in an analogous way as the one 
described above. In this simulation, the orthorhombic box was directly filled with water 
molecules. Moreover, this system was not heated at 473 K, since water molecules are 
much more mobile than sugars. The first two and four ns were not considered to compute 
the structural and dynamical properties presented in this paper for the 0-50 and 60 wt % 
systems, respectively. Table 1 summarizes some simulation data for the different systems 
considered in the present study. 

Table 1: System compositions (where N^, N5 and denote the number of lysozyme, 
sugar and water molecules, respectively), densities, and equilibration/simulation times for 
the different sugar concentrations (on a protein-free basis). Data corresponding to = 
wt % result from only one simulation of the lysozyme/pure water solution. 



(Wt %) 


N L /N 5 /N W 


density (g.< 


cm 3 ) 


Eq./Sim. time (ns) 


T M 


S 


T, M, S 





1/0/3800 


1.04 1.04 


1.04 


2/8 


37 


1/85/2800 


1.16 1.16 


1.15 


2/8 


50 


1/125/2400 


1.20 1.21 


1.20 


2/10 


60 


1/165/2100 


1.24 1.25 


1.24 


4/13 



3 Structural properties. 

3.1 Lysozyme 

3.1.1 Protein conformation 

The influence of sugars on the conformation of lysozyme has first been checked by comput- 
ing the root- mean-square deviation (RMSD) from the crystallographic structure [48], of 
either the C a carbon atoms of the protein backbone or of all protein atoms. The averaged 
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values for all studied systems are shown in table 2. The values for the pure water system 
are in fair agreement with those reported in previous studies [12,49], which employed the 
TIP3P water model [52]. The presence of sugars at the studied concentrations reduces 
conformational changes of lysozyme, as seen from the significant decrease of the C a and 
the all-atom RMSDs, which tends to be more important when disaccharide concentra- 
tion increases. A comparable effect has been observed by Lins et al. [12] for lysozyme 
in presence of trehalose at a concentration of ~ 18 wt % compared to the pure water 
system, where the all-atom RMSD of lysozyme in presence of trehalose was found about 
0.2 A lower. This reduction of the RMSDs probably arise from the slowing-down of the 
solvent dynamics observed at these concentrations (see section 4). No major difference is 
observed among the three sugars, whatever the concentration considered. This suggests 
that they have a similar influence on lysozyme from a structural point of view. This is 
expected in the framework of the preferential hydration hypothesis [22], which proposes 
that sugars are preferentially excluded from the protein surface. The conformation of 
lysozyme in the different solutions has also been characterized by means of its radius 
of gyration R g and its total solvent accessible surface area (SASA) (calculated with the 
program DSSP [53]). These parameters are given in the remaining part of table 2. The 
R g and SASA of lysozyme seem to be slightly larger in the presence of sugars. This may 
stem from the non-negligible interaction of surface residues of lysozyme with sugars, as 
will be shown in the section 3.2.1. Again, a clear distinction among the studied sugars 
does not emerge. These three parameters show that the structure of lysozyme in presence 
of sugars remains relatively close to that in pure water. 

Table 2: Parameters describing the structure of lysozyme : (i) Root-mean-square deviation 
(RMSD) from the crystallographic structure of C a carbon atoms or of all atoms, (ii) radius 
of gyration R g of lysozyme, and (iii) total solvent accessible surface area (SASA). Standard 
deviations from mean values are given in parentheses. 



(Wt %) 






RMSD (A) 


MA) 


SASA (A 2 ) 




Sugar 


Backbone 


All atoms 











1.40 


(0.13) 


2.27 (0.11) 


14.07 (0.07) 


7280 (144) 




T 


1.01 


(0.10) 


1.80 (0.08) 


14.21 (0.10) 


7380 (208) 


37 


M 


1.09 


(0.07) 


1.93 (0.06) 


14.20 (0.10) 


7310 (135) 




S 


1.06 


(0.06) 


1.84 (0.06) 


14.25 (0.05) 


7335 (122) 




T 


0.96 


(0.07) 


1.73 (0.06) 


14.14 (0.07) 


7167 (102) 


50 


M 


1.02 


(0.06) 


1.71 (0.06) 


14.21 (0.05) 


7290 (91) 




S 


0.97 


(0.11) 


1.66 (0.12) 


14.19 (0.07) 


7358 (189) 




T 


0.97 


(0.07) 


1.70 (0.07) 


14.26 (0.06) 


7437 (96) 


60 


M 


0.96 


(0.07) 


1.66 (0.08) 


14.28 (0.04) 


7325 (91) 




S 


0.94 


(0.06) 


1.63 (0.05) 


14.23 (0.05) 


7283 (73) 
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3.1.2 Protein fluctuations 



The influence of sugars on the internal motions of lysozyme has been investigated by 
the computation of atomic mean-square fluctuations (MSFs). Similarly to Maragliano et 
al. [54], MSFs were calculated by averaging over 250-ps time blocks to limit the effect of 
conformational changes during simulations. The figure 2 presents the MSFs of lysozyme 
main-chain atoms (C a , C, and N) for the different studied systems, averaged by residue 
after the removal of the overall translational and rotational motions of the protein using a 
least-square fitting procedure. These MSFs are in qualitative agreement with the exper- 
imental temperature factors of lysozyme in its tetragonal form [55]. The high diversity 
of fluctuations reflects the heterogeneity of local environments experienced by lysozyme 
residues. This is consistent with the three kinds of residues identified by Lumry [56] , based 
on temperature factors, and also with the knot /matrix classification in the description of 
Gregory et al. [38] from exchange rate distributions : residues with low and intermediate 
B-factors are identified to knot and matrix residues, respectively, whereas residues with 
the highest B-factors are located on the surface of lysozyme. In figure 2, five main parts 
of lysozyme (labelled from I to V) have rather high MSFs. These residues correspond 
to loops, 3-10 a-helices, and some residues of the /9-sheet. They are all located at the 
protein surface, and are therefore likely to form HBs with the solvent. Two zones of the 
lysozyme residue sequence (labelled A and B) have lower MSFs than the average. They 
correspond to an a-helix and a loop located in the core of lysozyme. These residues are 
not extensively accessible to the solvent (see section 3.1.3) and their motions are sterically 
constrained by the presence of other protein residues. The addition of sugars significantly 
reduces lysozyme fluctuations, but their distributions remain qualitatively the same for all 
sugar concentrations. The MSFs do not reveal clear differences among the three sugars, 
even if maltose tends to globally reduce slightly more protein fluctuations than trehalose 
and sucrose. This result appears in fair agreement with the figure la of ref. [57], where the 
heat capacity at constant pressure C p of the ternary lysozyme/maltose/water solution at 
40 wt % is found lower than the C p of the corresponding sucrose and trehalose solutions 
and may thus indicate that lysozyme is less flexible in presence of maltose than in presence 
of sucrose or trehalose. The strong decrease of lysozyme atomic fluctuations with addition 
of sugars suggests that its dynamics is slaved to that of the solvent (see section 4) i.e. 
governed by the solvent viscosity. This observation at 300 K implies that sugars may be 
able to shift the denaturation temperature T m of proteins towards higher temperatures, 
by reducing the amplitudes of local motions that might lead to denaturation. This is 
actually what has been observed on lysozyme and other proteins [20,57]. 
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Figure 2: Residue- averaged mean square fluctuations (MSFs) of atoms from the backbone 
of lysozyme (C a ,C,N) in the lysozyme/pure water solution and in the different ternary 
solutions (the overall translational and rotational motions of the protein have been re- 
moved using a least-square fitting procedure). The thick and thin black horizontal lines 
that appears in the top of the figure (in the wt % plot) denote a- and 3-10 a helices, 
respectively, whereas the thick grey lines indicate the /9-sheet. 



3.1.3 Fractional solvent accessibilies f sa 

The fluctuations of lysozyme residues clearly reveal distinct local environments, which 
may arise from differences in local packing and interactions with the solvent. Indeed, 
residues from the core should be more densely packed and less accessible to the solvent 
than residues from the surface, whose motions are also less sterically constrained. The 
time-averaged fractional solvent accessibilities [58] f sa of lysozyme residues have been cal- 
culated to know how they may interact with the solvent. The DSSP program has been 
used for this purpose [53]. f sa is equal to the ratio A prote i n / 'Af ree , where A prot ein is the 
SASA of a particular residue in the presence of the other surrounding residues of lysozyme, 
while Af ree is the related surface area of the free residue, i.e. without the presence of 
other amino acids. Figure 3 shows the f sa of the different lysozyme residues in pure water 
and in the 50 wt % trehalose solution (the distributions for the remaining solutions are 
similar and are thus not shown). These distributions are close to each other and further 
evidence that the presence of sugars do not modify significantly the native conformation 
of lysozyme. This is well in line with the results from Raman scattering measurements in 
the amide I region of lysozyme (1550-1750 cm -1 , see figure 5 of ref. [57]), which showed 
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that spectra in presence or in absence of sugars at 40 wt % are very similar to each other. 
This is expected since sugars are well-known protein stabilizers. Indeed, according to 
the preferential hydration hypothesis [22], sugars, and more generally osmolytes, desta- 
bilize much more the denatured state of proteins than their native state, and therefore 
lead to their stabilization against the denaturation process. A more important point is 
that the distributions mimic those of MSFs (see figure 2), and thus confirm that residues 
with low solvent accessibilities fluctuate much less than residues exposed to the solvent. 
Figure 3c presents the distributions probabilities of f sa . These distributions show a bi- 
modal behavior : (i) a first relatively sharp contribution for f sa lower than about 0.15, 
and (ii) a second much broader band, centered at ~ 0.3. These two contributions have 
been attributed to lysozyme core and surface residues, respectively. The sharper contri- 
bution from core residues may indicate that the environments they experience is much 
homogeneous than that of surface residues. This again shows that surface residues may 
exhibit larger motions than core residues, because sterical constraints imposed by other 
protein residues are lower. This result also suggests that the exposition of residues to the 
solvent may have important consequences on their dynamical properties (see section 4). 
A detailed analysis of the solvent properties is therefore needed to understand how sugars 
may preserve proteins against denaturation. 
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Figure 3: Fractional solvent accessible surface areas f sa as a function of the residue number 
sequence for the lysozyme/pure water (a) and for the 50 wt % trehalose (b) solutions 
(smoothed curves computed with the Loess algorithm [59] serve as guides to the eye). 
The related probability distribution functions are displayed in (c). 
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3.2 Protein-solvent interactions. 
3.2.1 Protein-solvent HBs 

The way sugars interact with biomolecules, and in particular with globular proteins and 
membrane phospholipid bilayers, has been largely debated over the last two decades [6]. 
According to the water replacement hypothesis, sugars form HBs with biomolecules by 
substituting to their hydration water and may then preserve their native structure. In 
the present study, the HBs formed at the lysozyme-solvent interface have been analyzed. 
Two molecules were considered to be H-bonded if the donor-acceptor distance do a is less 
than 3.4 A and if the D-H- ■ -A angle is larger than 120 deg. [60]. The table 3 summarizes 
the time-averaged numbers of HBs formed between lysozyme and solvent molecules. In 
the investigated concentration range, sugars form between about 10 and 20 % of the total 
number of protein-solvent HBs, as indicated in table 3. This is much fewer than the 
number of HBs they could form, assuming that water and sugar molecules were equally 
able to form HBs with lysozyme. Indeed, assuming that water is able to form 4 HBs (2 as 
donor and 2 as acceptor) and that each hydroxyl group of the disaccharides may optimally 
form 3 HBs (1 as donor and 2 as acceptor), the proportion of sugar-lysozyme HBs among 
solvent-lysozyme HBs should be about 15.4, 23.8 and 32.0 % for the 37, 50, and 60 wt % 
solutions, when sugar-lysozyme HBs via ring and glycosidic oxygen atoms of sugars are 
neglected. Therefore, sugars are preferentially excluded from the surface of lysozyme, as 
suggested by the preferential hydration hypothesis [19,22]. In addition, the total number 
of protein-solvent HBs tends to slightly decrease when sugar concentration increases. This 
suggests that sugars are not able to substitute perfectly to water molecules. They indeed 
cannot access as easily as water molecules to polar groups of lysozyme because of sterical 
hindrance and topological constraints. 

The number of disaccharide-protein HBs and the number of water molecules shared 
between lysozyme and disaccharides increase with sugar concentration, since sugar-protein 
interactions are more likely. Slight differences are observed among the three sugars. First, 
sucrose molecules seem to form a lower number of HBs with lysozyme in the 50 wt % 
solution. This difference may partly arise from the fact that the HB network of sucrose 
is less developed than those of trehalose and maltose, which have almost percolated at 
this concentration (see figure 5d). It is also likely to occur from a lack of configurational 
sampling of lysozyme-sugar interactions. Indeed, sucrose form more comparable numbers 
of HBs with lysozyme at lower or higher concentrations. Consequently, care must be taken 
when comparing the results for the three sugars. Moreover, it does not seem that trehalose 
forms more HBs with lysozyme than sucrose and maltose at the studied concentrations, 
although this may be the case in the completely dehydrated case. This result would 
be consistent with the simulation results of Lins et al. [12] which showed that trehalose 
clusters at the protein surface and does not expel the water molecules closest to the protein 
surface. It would therefore confirm the preferential hydration hypothesis proposed by 
Timasheff et al. [21,22], which suggests that trehalose is more excluded - in the relatively 
diluted solutions - than other osmolytes, and would thus increase the thermodynamical 
stabilization of the proteins compact native state relative to their extended denaturated 



12 



state [3]. The larger hydration number of trehalose (see figure 5a) would then explain why 
it does not interact more than do maltose and sucrose with lysozyme in our simulations. 
These results also appear in line with those of a recent simulation study [36] , which shows 
that maltose forms a larger number of HBs with a phospholipid membrane than trehalose, 
at a 2 m concentration, at both 325 and 475 K. 

Finally, the water replacement hypothesis does not seem confirmed by these results, 
since sugar-protein HBs are not able to substitute perfectly to water-protein HBs, prob- 
ably because of topological constraints and excluded volume effects. Nevertheless, sugar- 
lysozyme HBs relax slower than water-lysozyme HBs as a consequence of the lower mo- 
bility of sugars, so that a fewer number of sugar-protein HBs may actually increase the 
stabilization of proteins. 

Table 3: Mean numbers of water-lysozyme and disaccharide- lysozyme HBs, n//s(W-L) 
and tih\b(D-L), respectively, proportion of disaccharide- lysozyme HBs among solvent- 
lysozyme HBs (% HB (D-L) = n HB (D-L)/(n HB {D-L)+n HB (W-L))) and number n w (L- 
W-D) of water molecules which are simultaneously hydrogen-bonded to lysozyme and 
disaccharides. Standard deviations from mean values are given in parentheses. 



(Wt %) 




n HB (W-L) 


iWD-L) 


% HB (D-L) 


n^(L-W-D) 







329.0 (5.1) 










T 


291.8 (6.0) 


33.9 (3.2) 


10.4 (1.0) 


46.2 (4.4) 


37 


M 


285.7 (5.6) 


45.4 (4.7) 


13.7 (1.4) 


54.2 (4.0) 




S 


277.1 (5.6) 


35.5 (3.3) 


11.4 (1.1) 


49.4 (3.6) 




T 


261.4 (4.5) 


50.9 (3.3) 


16.3 (0.9) 


62.4 (3.8) 


50 


M 


265.1 (4.2) 


54.9 (4.1) 


17.2 (1.2) 


62.2 (6.7) 




S 


281.5 (5.8) 


34.1 (3.7) 


10.8 (1.1) 


57.3 (4.3) 




T 


259.7 (5.6) 


57.0 (4.2) 


18.0 (1.3) 


73.0 (3.2) 


60 


M 


247.2 (6.1) 


62.1 (4.1) 


20.1 (1.4) 


74.4 (5.1) 




S 


243.6 (6.6) 


64.0 (4.1) 


20.8 (1.4) 


65.7 (3.9) 
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3.2.2 Preferential hydration 



Solvent molecules do not form HBs with apolar groups of lysozyme. Therefore, the number 
of protein-solvent HBs does not describe exhaustively protein-solvent interactions. We 
have thus characterized the relative local distribution of water molecules around lysozyme 
in a similar way to Cottone et al. [13]. We have indeed computed the time-averaged 
normalized ratio gN,Ow=now /( n ow + n s) (?) / (N ow / (N ow + N os )), where n w and 
n os are the local numbers of water oxygen atoms and sugar hydroxyl group oxygens, 
respectively, located at a minimum distance r from any heavy atom of lysozyme, and 
Now and Nos denote the total numbers of water oxygen atoms and sugar hydroxyl group 
oxygen atoms, respectively, in the simulation box. This ratio is greater than one in the 
close proximity of the protein surface if a given sugar is preferentially excluded from the 
protein surface (in other words, if the protein is preferentially hydrated). Conversely, the 
sugar preferentially interacts with the protein if this ratio is lower than one. This ratio 
is represented in the figure 4 for the different ternary studied systems. Lysozyme clearly 
appears more and more preferentially hydrated when increasing sugar concentration. At 
distances larger than ~ 5 A there is a slight water depletion which results from the 
presence of sugars. This preferential hydration probably arise from excluded volume 
effects and is consistent with the relatively low proportion of sugar-lysozyme HBs relative 
to water- lysozyme HBs (see table 3). This may confirm that sugars are good protein 
stabilizers in the framework of the preferential hydration hypothesis [22]. If we exclude 
the sucrose solution at 50 wt %, for which there may be shortcomings, we note that 
trehalose preferentially hydrates lysozyme slightly more than do maltose and sucrose. 
This seems to be especially true for the hydration of apolar groups of lysozyme, which 
is apparent in the small peak located around 3.7 A . Trehalose hydroxyl oxygens were 
indeed systematically more excluded from these groups than those of maltose and sucrose 
(data not shown). This may arise from the larger hydration number of trehalose (see 
figure 5a), which would prevent it from remaining close to apolar groups of lysozyme. 
Cottone et al. [13] have also shown that sucrose and trehalose at a concentration of 89 
wt % preferentially hydrate carboxy-myoglobin (MbCO), but sucrose to a larger extent 
than trehalose. It is possible that the effect of trehalose is somewhat different at this 
much more elevated concentration, as exemplified by the differences between the water 
replacement [6] and the preferential hydration [22] hypotheses. 
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Figure 4: Normalized fraction of water oxygen atoms QN,Ow=now / ( n ow + 
no s) I '(Now I '(Now + Nos)) as a function of the minimal distance to any heavy atom 
from lysozyme. now and nos denote the local numbers of water and sugar hydroxyl 
oxygen atoms, respectively. Now an d Nos correspond to the total numbers of water and 
sugar hydroxyl oxygen atoms found in the simulation box, respectively. 



3.3 Solvent structure. 

An important issue of the bioprotection phenomenon involves the role of sugar-water in- 
teractions. Indeed, the stabilizing effect of solutes is sometimes thougth to depend directly 
on their effect on the HB network of water, although this is controversed [61]. Therefore, 
the larger hydration number of trehalose compared to maltose and sucrose [23, 24, 39] may 
explain its enhanced stabilizing effect. In the present solutions, the interaction of sugars 
with water molecules is rather analogous to that of binary water/sugar solutions [39]. 
First, the figure 5a reveals that the hydration numbers of sugars are lower than in our 
previous simulation study of binary sugar/water mixtures [39], even if we exclude the 
sugars hydrogen bonded to lysozyme (data not shown). The analysis of additional sim- 
ulations of dilute binary sugar/water solutions at a concentration of 4 wt % (1 sugar 
with 512 water molecules) using the CSFF carbohydrate force field [42] suggests that the 
observed differences between binary and ternary solutions do not arise from the different 
force fields, Ha et al. [62] and CSFF [42], used for representing sugars in ref. [39] and in 
the present study, respectively. It could thus rather be attributed to the preferential ex- 
clusion of sugars, which raises their concentrations within their corresponding accessible 
volumes, thus reducing their hydration numbers. Trehalose is found to have a slightly 
larger hydration number than maltose and sucrose, with the exception of the 60 wt % so- 
lution, where it is found slightly less hydrated than maltose, as shown in the figure 5a. It 
is possible that the hydration behavior of sugars in ternary solutions cannot be perfectly 
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extrapolated from that in binary aqueous solutions, because there should exist for sugars 
at high concentrations a compromise between having a large hydration number and being 
highly preferentially excluded from the protein surface. Secondly, the numbers of sugar 
intramolecular HBs, shown in figure 5b, are also in qualitative agreement with our previ- 
ous study : trehalose is found to form a lower number of intramolecular HBs than does 
sucrose, maltose being intermediate, whatever the concentration considered. This would 
explain why the hydration number of trehalose is slightly larger than those of maltose 
and sucrose, since hydroxyl groups involved in these HBs no longer remain available for 
interacting with water molecules. These numbers definitely depend on the force field used 
and are systematically found lower in the present study. This is consistent with the larger 
radii of gyration of sugars with the CSFF force field [42] (3.49 A, 3.52 A, 3.41 A for T, 
M, S, respectively) than with the Ha et al. force field [62] (3.40 A, 3.46 A, 3.30 A for 
T, M, S, respectively) in the dilute binary solutions (4 wt %). This result suggests that 
sugar conformations are sligthly more extended and that their rings are slightly further 
from each other with the CSFF force field. Moreover, the SASAs of sugars computed 
according to the Lee and Richards method [63] with a probe radius of 1.4 A are shown 
in figure 5c. The SASA of trehalose is found larger than that of sucrose, that of maltose 
being intermediate. This property may explain why trehalose interacts more with water 
than do maltose and sucrose [23,24]. Besides, the presence of a furanose ring in sucrose 
(see figure 1) probably explain why several of its properties differ from those of maltose 
and trehalose. 

Finally, the percolation of the HB network of sugars was suggested to play a major 
role in the bioprotection phenomenon, and differences between the three sugars have 
been observed [39]. This is particularly true for sugar concentrations above 50 wt %, 
where sugar-sugar interactions become important. The normalized mean sugar cluster size 
< ris > /Ns has been computed in the same way as in our previous work [39] to know if the 
percolation of the sugar HB network occurs in the same concentration range in presence of 
the protein. These sizes are presented in the figure 5d. A transition clearly appears in the 
[37-50] wt % concentration interval and the HB network of sugars has almost completely 
percolated at 60 wt %, in good agreement with our previous results [39]. Sucrose on 
one hand and maltose and trehalose on the other hand show quite different behaviors. 
Indeed, the percolation transition in sucrose systems is shifted to higher concentrations 
compared to maltose and trehalose ones. In other words, a higher concentration of sucrose 
is needed to reach a given cluster size < ns >■ This is well in line with the higher number 
of intramolecular HBs in sucrose (see figure 5b), as well as its slightly lower accessible 
surface area due to the presence of a furanose ring (see figure 5c). This behavior may 
then explain why the dynamics of lysozyme is less modified by the presence of sucrose 
than by the presence of maltose or trehalose, as shown in the section 4. Furthermore, 
trehalose and maltose, which are topologically closer to each other than to sucrose (see 
figure 1), behave in a similar way. Nonetheless, it seems that the HB network of maltose 
tends to percolate at a lower sugar concentration than that of trehalose, as observed 
previously [39]. This may originate both from the less symmetric conformation of maltose 
and from its lower affinity for water molecules. In our previous study, maltose formed 
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systematically larger clusters than trehalose, whereas this not always remains true in 
the present study. This difference may originate from the presence of the protein, with 
which maltose seems to interact more strongly than trehalose does and thus less with 
other maltose molecules. Interestingly, the mean sugar cluster sizes < ns > /Ns seem 
to be consistent with the solubility of the three disaccharides. Indeed, the experimental 
solubility at room temperature of sucrose is significantly larger than that of trehalose, 
whereas the solubility of maltose is slightly lower than that of trehalose (see fig. 4 of 
ref [64]). Recently, Giuffrida et al. [65] have suggested that this enhanced tendency of 
maltose molecules to form clusters could be related to their larger dipole momentum. 
Indeed, from ab initio calculations at the BLYP/6-31G** level, they obtained values 
of 5.2, 2.5 and 1.5 Debyes for the crystallographic structures of maltose, sucrose, and 
trehalose, respectively. 




40 45 50 55 60 

Disaccharide concentration (wt %) 



Figure 5: Hydration number n H (a), normalized mean number of intramolecular HBs 
< n HB > intra /Ng (b), solvent accessible surface area (SASA) (c) and normalized mean 
sugar cluster size < ns > /Ns (d) of trehalose (solid line), maltose (dashed line), and 
sucrose (dotted line) as a function of their concentration. 



4 Dynamical properties. 

It is well known that proteins exhibit extremely diverse motions (from local atomic to 
subunit motions) which occur on time scales that span many orders of magnitude [66, 67]. 
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The broad distribution of the MSFs and of the f sa of lysozyme residues actually suggest 
that their structural relaxations are very heterogeneous, and thus that lysozyme dynamics 
is complex. A recent study of hydrated lysozyme (h = 0.43) indeed shows that the 
distribution of the effective rotational correlation times of methyl groups is very broad [68]. 
Furthermore, the effect of sugars on the MSFs of lysozyme suggest that they slow down 
its diffusive motions. This would imply that they make lysozyme less sensitive to high 
temperatures by preventing large-scale motions that might lead to denaturation. In order 
to get a deeper insight into the influence of sugars on the relaxational behavior of lysozyme, 
we have computed incoherent intermediate scattering functions Si nc (Q,t) for each residue 
of lysozyme. Si nc (Q, t) is defined as : 

^(Q,^(i:c et4[Mt) " Mo)i } 

a 

where b a ^ nc and r denote the incoherent scattering length and the vector position of a 
given atom a, respectively, and Q is the wavevector - the brackets mean averaging over 
every time origin of simulations. These functions represent the time Fourier transform 
of dynamic structure factors Si nc (Q,cu), obtained in neutron scattering experiments. We 
have chosen to probe the dynamics of lysozyme at a wavevector of 2.29 A -1 , which 
corresponds to the position of the first peak in the static structure factor So-o of pure 
water and which was used in our study of binary sugar/water solutions [50]. Charac- 
teristic relaxation times r may then be defined as the decay times from 1 to 1/e of the 
different S inc (Q,t). 

The figure 6 shows the relaxation frequencies 1/r of residues of lysozyme in pure 
water and in the different trehalose solutions. These frequencies clearly appear diverse 
and, similarly to MSFs, they seem to follow the fractional solvent accessibilies f sa of 
residues of lysozyme (see figure 3), since residues which are more exposed to the solvent 
relax faster than those which are buried in the protein core. 

The addition of sugars at the studied concentrations induces a strong reduction of 
the relaxation frequency of lysozyme residues. This slowing-down of the dynamics of 
lysozyme evidences that the higher the solvent viscosity, the lower the diffusive (or a-like) 
motions of lysozyme, in the studied concentration range. This would confirm the basic 
assumption that the structural or a relaxation of the solvent - which is related to its 
viscosity - determines the influence of the solvent on the dynamics of proteins [69]. 

The dynamical slowing-down induced by sugars is more visible in the figure 7a, which 
shows the probability distribution functions of the relaxation times r of lysozyme residues 
in the pure water and in trehalose solutions (the distribution functions for maltose and 
sucrose solutions exhibit comparable behaviors and are not shown for clarity reasons). A 
significant shift towards larger relaxation times as well as a broadening of distributions are 
observed when the sugar concentration increases. This may reflect the broadening of the 
energy barriers experienced by lysozyme residues resulting from an increased dynamical 
coupling between the solvent and the protein. Furthermore, the correlation between f sa 
and relaxation times r remains relatively well established at the highest studied concen- 
tration (60 wt %), as shown in figure 7b, where is represented the mean relaxation time 
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Figure 6: Relaxation frequencies 1/r of the residues of lysozyme in pure water and in the 
different trehalose solutions (smoothed curves computed with the Loess algorithm 
serve as guides to the eye). 



of lysozyme residues as a function of their fractional accessible surface area f sa , for f sa 
up to of about 0.6 (given the low number of residues with a f sa larger than 0.6 as shown 
in figure 3, the related statistical errors are very large). This figure shows that sugars 
influence both the dynamics of core and surface residues. It seems therefore possible to 
modulate the internal dynamics of proteins by changing the solvent viscosity with the 
addition of sugars. 

Figure 7b also shows that the relative influence of sugars on the dynamics of lysozyme 
does not change too much over the entire f sa range. Therefore, a straigthforward com- 
parison between sugars is possible by considering the mean relaxation time of the whole 
lysozyme. The mean relaxation times of sugar and water molecules have also been com- 
puted to interpret the effect of sugars on lysozyme and solvent dynamics. The figure 8 
shows that the dependences on the disaccharide concentration of the mean relaxation 
times of lysozyme, sugars and water molecules exhibit a similar slope change at a concen- 
tration close to 50 w %. This may be interpreted when considering the percolations of 
the sugar HB networks, which occur at a concentration of about 40-50 wt % in both the 
sugar binary [39] and ternary solutions (see figure 5d). At these relatively high concen- 
trations, the hydration layers of sugars do not surround completely disaccharides because 
of sugar-sugar HBs, which strongly influence their dynamics. These dynamical changes 
have been observed experimentally by Rampp et al. [70] for a series of carbohydrates - 
including sucrose and trehalose - at concentrations above 50 wt %. They also appear 
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Figure 7: (a) Probability distribution functions of the relaxation times r of lysozyme 
residues in the pure water /lysozyme solution and in the different trehalose solutions, (b) 
Averaged relaxation times of lysozyme residues r as a function of their fractional solvent 
accessibilities f sa for the pure water solution and for the 60 wt % trehalose, maltose and 
sucrose solutions. 

in the concentration dependences of the diffusion coefficient of sucrose and trehalose in 
aqueous solutions obtained by Ekdawi-Sever et al. by NMR measurements [71]. Sugar- 
sugar interactions are not large enough at 37 wt % for detecting significant differences 
between sugar solutions. But at 50 wt %, the HB network of maltose and trehalose is more 
developed than that of sucrose, so that each species relaxes more slowly in the maltose 
and the trehalose solutions than in the sucrose one. Finally, the HB network of sugars 
has almost percolated in the 60 wt % solutions. Their dynamics are then dominated by 
the relaxation of the HB network of sugars and differences among sugars emerge. Sucrose 
clearly relaxes faster than trehalose and maltose, due to its lower ability to form clusters. 
This is well in line with experimental studies, which show that the dynamics of sucrose is 
faster than that of maltose and trehalose [23]. Moreover, the relaxation times of maltose 
seem to be somewhat longer than those of trehalose. This appears consistent with the 
lower heat capacity at constant pressure C p of the ternary lysozyme/maltose/water solu- 
tion at 40 wt % in comparison with the corresponding sucrose and trehalose solutions [57], 
which suggests that the ternary maltose solution is more viscous. This result also seems 
to confirm the peculiar behavior of maltose observed by Giuffrida et al. [65], where the 
matrix dynamics was found the most reduced in the dry maltose sample. Experimen- 
tal measurements suggest that binary trehalose/ water solutions are sligthly more viscous 
than maltose ones [23,72]. Therefore, the present results may underline the influence of 
the protein on the dynamics of the solvent (see e.g. ref. [14] and ref. [73]), which may 
arise from (i) the roughness of the protein surface, (ii) a decrease of the dimensionality 
of the solvent diffusion at the protein-solvent interface and (iii) strong solvent-protein in- 
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teractions with polar and charged groups of the protein. Given that the conformation of 
lysozyme is similar in trehalose and maltose solutions, the stronger interaction of maltose 
with lysozyme (see table 3) may imply that its dynamics is more slowed down than that 
of trehalose. 

Furthermore, a direct comparison between the expected biopreservative efficiencies of 
the three studied carbohydrates in the light of the present results is not straightforward 
since maltose is a reducing sugar, unlike trehalose and sucrose. It may thus react with 
lysine and arginine residues of lysozyme and other proteins to form carbohydrate adducts, 
especially at high temperatures, via a complex browning pathway called Maillard reac- 
tion [3] (which cannot occur in MD simulations). This may explain why maltose is rarely 
used experimentally for protein stabilization, in contrast to trehalose and sucrose, which 
are natural bioprotectants. Maltose was actually found to induce the smallest shift of the 
denaturation temperature T m of lysozyme at a concentration of 40 wt %, as measured in 
modulated differential scanning calorimetry and Raman scattering investigations [57]. 

In addition, figure 8 reveals that the protein, sugar and water dynamics are slower in 
the trehalose solutions than in the sucrose solutions for concentrations above 40-50 wt %. 
These results are in agreement with the tighter coupling of carboxy-myoglobin (MbCO) 
with trehalose-water matrices than with sucrose-water matrices observed experimentally 
by Giuffrida et al. by means of Fourier transform Infrared (FTIR) measurements [74]. 
This enhanced slowing-down presumably arises from the larger sugar-sugar interactions in 
the trehalose solutions at these concentrations, as seen in figure 5d. Moreover, the larger 
hydration number of trehalose compared to sucrose allows a more important dynamical 
coupling with water molecules. Therefore, trehalose would have a better preservation 
efficiency because of the quasi-absence of intramolecular HBs in its solution conforma- 
tion (see figure 5b), which allows it to interact more strongly with both water and sugar 
molecules than sucrose does. On the other hand, this would explain its higher solubility 
than (true) maltose [64] . This suggests that the greater homogeneity [39] of intermolecu- 
lar interactions in the trehalose solutions (water-trehalose and trehalose-trehalose) might 
lead to their better preservation efficiency. This may explain why trehalose was found 
to induce the largest shift of T m and to increase the most the temperature interval of 
the denaturation process AT of lysozyme [57]. The shift of T m was indeed found to be 
of 9.4±0.3 K, 4.2±0.3 K and 2.9±0.3 K in the trehalose, sucrose and maltose solutions, 
respectively, by Raman scattering. Similarly, the temperature domain AT of the denat- 
uration process was extended by 0.9±0.2 K, 0.4±0.2 K and 0.5±0.2 K in the trehalose, 
sucrose and maltose solutions, respectively (see table 2 of ref. [57]). 

Finally, it must be pointed out that the slowing-down of lysozyme dynamics follows 
the slowing-down of the solvent induced by the percolation of the sugar HB network, 
but not necessarily the number of lysozyme- sugar HBs. Indeed, sucrose forms at 60 wt 
% a larger number of HBs with lysozyme than trehalose, but lysozyme relaxes faster 
in the sucrose solution. This would imply that the slowing-down of the solvent plays 
a more important role on the dynamics of lysozyme than specific sugar-lysozyme HB 
interactions in the studied concentration range. If we only compare trehalose and sucrose, 
this suggests that the hypothesis of Green and Angell [27], as well as the preferential 
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hydration hypothesis [19, 22] might be valid for this kind of solutions. In presence of 
sugars, the unfolding process would require a reorganization of the HB network of sugars, 
easier in the sucrose solutions than in the trehalose ones. 
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Figure 8: Mean relaxation times r of lysozyme (top), disaccharides (middle) and water 
(bottom) as a function of the disaccharide concentration. 
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5 Conclusions 



This article reports the results from molecular dynamics simulations of lysozyme/disaccha- 
ride/water solutions for intermediate concentrations of disaccharides (37-60 wt % on a 
protein- free basis). This concentration range was chosen because it corresponds to the 
range of concentrations where the HB network of sugars was shown to percolate in our 
previous study of disaccharide/water solutions [39]. Several parameters of lysozyme like 
its root-mean-square deviations (RMSD) from the crystallographic structure, its radius 
of gyration and its SASA indicate that its conformation in presence of sugars remains 
very similar to that in pure water. Nevertheless, structural fluctuations of lysozyme as 
seen from mean-square fluctuations (MSFs) are much reduced in presence of sugars. This 
reduction was shown to be quite homogeneous throughout the protein and follows the 
fractional solvent accessibilities f sa of lysozyme residues. This points out the influence 
of the solvent on the internal dynamics of lysozyme. The analysis of the interactions of 
sugars with lysozyme via HBs as well as the relative concentration of water oxygen atoms 
around lysozyme suggest that lysozyme is preferentially hydrated, well in line with other 
studies [12, 13]. When comparing the three sugars, it seems that trehalose is slightly more 
excluded than maltose and sucrose, although the observed differences may be not signifi- 
cant. This could be interpreted by the larger hydration number of trehalose compared to 
that of maltose and sucrose [23,24]. Furthermore, the number of intramolecular HBs of 
trehalose was found lower than that of maltose and sucrose. Interestingly, trehalose was 
also found to have the largest SASA. This may explain why it interacts more with water 
than maltose and sucrose. In addition, the HB network of sugars was shown to percolate 
in the studied concentration range. Sucrose forms smaller clusters than trehalose and 
maltose, probably because of its larger number of intramolecular HBs. Dynamical prop- 
erties of lysozyme have been analyzed via incoherent intermediate scattering functions 
Sinc{Q,t). The relaxation times of lysozyme residues mimic both MSFs and f sa distri- 
butions and a broadening of the distributions of relaxation times of lysozyme residues is 
observed when increasing sugar concentration. This could reveal an enhanced dynamical 
coupling between solvent and lysozyme. Moreover, a similar increase of lysozyme, sugar 
and water relaxation times is observed when increasing the sugar concentration from 37 
to 60 wt % and may arise from the percolation of the HB network of sugars, as suggested 
from our previous results on binary sugar/water solutions [39]. Since sucrose forms smaller 
clusters than trehalose and maltose, it appears to slow down the dynamics of lysozyme 
less strongly at concentrations above ~ 40-50 wt %. This appears consistent with experi- 
mental results, which sugggest a looser protein-solvent coupling in sucrose-water matrices 
than in trehalose-water matrices of various water contents [74]. In contrast, maltose slows 
down more efficiently the dynamics of lysozyme because it forms larger clusters and in- 
teracts more strongly with lysozyme than trehalose. This seems consistent with the lower 
heat capacity at constant pressure C p of lysozyme/maltose/ water solutions at 40 wt % 
compared to that of trehalose and sucrose solutions [57] and with the peculiar behavior of 
maltose observed experimentally by Giuffrida et al. [65]. 

Our simulations only cover a limited concentration range of disaccharides at ambient 
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temperature (300 K), so that lysozyme is not submitted to temperature and/or desiccation 
stresses. Furthermore, the length of simulations is too short to account for phenonema 
which occur on much longer time scales in real systems, such as vitrification, ice formation 
or sugar crystalline phase transitions. Nevertheless, our results suggest that the dynam- 
ical slowing-down of lysozyme induced by the presence of sugars may stabilize globular 
proteins against thermal denaturation. In addition, our results indicate that lysozyme 
remains preferentially hydrated in the studied concentration range, since the substitution 
of water-lysozyme HBs by sugar-lysozyme HBs is rather limited. 

The comparison of the expected biopreservative behaviors of the three studied disac- 
charides on lysozyme is not straightforward since maltose is a reducing sugar and may 
thus alter the structure of lysozyme when chemically reacting with its lysine and arginine 
residues via the Maillard reaction. Consequently, the comparison of sucrose and trehalose, 
which are naturally synthesized osmolytes, is more relevant in the context of preserving 
globular proteins. The results of this study indicate that trehalose has a stronger influ- 
ence on the dynamics of water and lysozyme than sucrose and may therefore preserve 
more efficiently the native structure of lysozyme at higher temperatures, particularly at 
concentrations above ~ 40-50 wt %. These differences appear to primarily stem from 
their different topologies (see figure 1), which cause the formation of intramolecular HBs 
in sucrose and - almost - none in trehalose. The main consequences are (i) a lower number 
of water molecules with which sucrose molecules may interact and (ii) the formation of 
sucrose clusters of lower sizes, which would have thus a reduced influence on the solvent 
and lysozyme dynamics. Even though these differences between trehalose and sucrose ap- 
pear relatively small at room temperature where lysozyme is in its compact native state, 
they may increase during the unfolding process of lysozyme and then play a major role. 
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